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Abstract. We show that Bose-Einstein condensates in optical lattices with broken 
time-reversal symmetry can support chiral edge modes originating from nontrivial 
bulk excitation band topology. To be specific, we analyze a Bose-Hubbard extension 
of the Kaldane model, which can be realized with recently developed techniques of 
manipulating honeycomb optical lattices. The topological properties of Bloch bands 
known for the noninteracting case are smoothly carried over to Bogoliubov excitation 
bands for the interacting case. We show that the parameter ranges that display 
topological bands enlarge with increasing the Kubbard interaction or the particle 
density. In the presence of sharp boundaries, chiral edge modes appear in the gap 
between topological excitation bands. We demonstrate that by coherently transferring 
a portion of a condensate into an edge mode, a density wave is formed along the edge 
owing to an interference with the background condensate. This offers a unique method 
of detecting an edge mode through a macroscopic quantum phenomenon. 
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1. Introduction 


Topological insulators and superconductors have attracted great attention in recent 
years for their rich variety of quantized responses and robust gapless edge states 
originating from nontrivial topology of bulk Bloch bands [0, |^, |^. A prototype of 
topological insulators is an integer quantum Hall system |^, which exhibits a quantized 
Hall conductivity a^y = —{e‘^/h)C, where C is the sum of the Chern numbers of 
occupied bands 0. The gapless edge states are characterized by \C\ sets of chiral 
modes propagating clockwise (counterclockwise) along the system’s edge for G > 0 
(C < 0) B 0. Here a nonzero value of C is caused by the breaking of time- 
reversal symmetry due to a magnetic field. The discovery of Z 2 topological insulators 
B.BE3O.0.0. P!4| has opened up a new avenue for realizing a topologically nontrivial 
structure in Bloch bands through spin-orbit coupling, without breaking time-reversal 
symmetry. Topological superconductors have been shown to exhibit exotic edge states 
consisting of Majorana fermions, which are protected by particle-hole symmetry [|T^, [l^. 
A unified understanding of these topological phases has been achieved with a topological 
periodic table, where such phases are systematically classihed for quadratic fermionic 
Hamiltonians in different dimensions and symmetry classes [^, 0. 

Ultracold atomic systems have recently emerged as a new platform for exploring the 
physics of topological phases, especially owing to ongoing experimental developments 
for engineering synthetic gauge fields |^, ^ which can be used to produce such 
states. Different schemes have been proposed and implemented to create nearly 
uniform magnetic helds in continuum |^, optical lattices |^, ^ and 

synthetic dimensions B00- Furthermore, the Haldane model in which non¬ 


uniform fluxes pierce through the system, has been realized using fermionic atoms in a 
periodically modulated honeycomb optical lattice (see Refs. ^ ^ ^ Q for 
early theoretical proposals for realizing the same or related models). The Haldane model 
is a prototypical example of Hamiltonians that exhibit topologically distinct regimes 
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characterized by the Chern numbers C± associated with upper (+) and lower (—) Bloch 
bands. The phase diagram of this model has been vindicated experimentally using 
momentum-resolved interband transitions |^. Another recent remarkable achievement 
has been an interferometric measurement of the vr Berry flux in the momentum space 


of a honeycomb lattice |36 


A notable feature of atomic systems is that one can study the effect of quantum 
statistics. For example, by implementing the technique of Ref. bosonic atoms, 

one can realize a bosonic counterpart of the Haldane model. In the noninteracting case, 
topological properties of Bloch bands do not depend on quantum statistics. For weakly 
interacting bosons in optical lattices, Bogoliubov excitation bands give the elementary 
excitations of Bose-Einstein condensates (BEC). It is then interesting to ask how the 
topological properties of Bloch bands are carried over to those of Bogoliubov excitation 
bands in the interacting case. 

Band topology of bosonic or classical vibrational modes has been studied previously 
in photonic |^, |^, |^, phononic ^ magnonic 


45 


47| , and 

polaritonic excitations. Nontrivial Chern numbers of bulk excitation bands give rise 
to in-gap chiral edge modes, as observed experimentally in photonic systems ||^, ^ . 

Ultracold bosonic atoms in optical lattices are expected to offer a unique platform for 
the studies of band topology because of the high controllability of such systems and a 
potential combination with the macroscopic quantum nature of BECs. 

In this paper, we study the topological properties of Bogoliubov excitation bands 
in BECs in optical lattices with broken time-reversal symmetry, using a Bose-Hubbard 
extension of the Haldane model (Haldane-Bose-Hubbard model). We show that the 
topological properties of the Bloch bands in the noninteracting case are smoothly 
carried over to those of the Bogoliubov excitation bands in the interacting case. 
Furthermore, the parameter ranges that exhibit nontrivial band topology enlarge with 
increasing the Hubbard interaction or the particle density (see Fig. In the presence of 
sharp boundaries, chiral edge modes appear in the gap between topologically nontrivial 
excitation bands. We demonstrate that by coherently transferring a portion of the 
condensate into an edge mode, a density wave is formed along the edge owing to 
an interference with the background condensate. This property can be used as a 
macroscopically enhanced signature of an edge mode. 

We note that Vasic et al. have recently studied the ground-state phase diagram 
of the Haldane-Bose-Hubbard model, predicting the emergence of uniform and chiral 
BEC phases and plaquette Mott insulators with loop currents. While the Bogoliubov 
excitation bands in the uniform and chiral BEC phases have also been studied, their 
topological properties such as Chern numbers and associated edge modes have not been 
analyzed in detail. Our work addresses such topological properties of excitations in 
a uniform BEC phase, clarihes the parameter ranges showing topologically nontrivial 
bands in the presence of interactions, and demonstrates a unique method of detecting 
an edge mode though a macroscopic quantum interference. We also note that strong 
correlation effects on the band topology have been studied in the spin-^ fermionic 
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Haldane-Hubbard model in Refs. |52 

The rest of the paper is organized as follows. In Sec. |^, we describe our model, and 
review the band structure in the noninteracting case. In Sec. |^, we present a Bogoliubov 
theory for homogeneous condensates with weak repulsive interactions, and analyze the 
topology of Bogoliubov excitation bands. In Secs. ^ and ||, we analyze the ground state 
and excitations of inhomogeneous condensates by using a Bogoliubov-de Gennes theory. 
After describing the basic formalism in Sec. we present numerical results for box and 
harmonic traps in Sec. |^. In Sec. P, we present a summary of this paper and discuss an 
outlook for future studies. 


2. Model and band structure 

In this section, we hrst describe our model—a Bose-Hubbard version of the Haldane 
model in a honeycomb lattice. We then review the band structure in the noninteracting 
case, and discuss the single-particle ground state into which bosons condense at zero 
temperature. 


2.1. Bose-Hubbard Hamiltonian 


We consider bosonic atoms in an optical lattice, which are well described in the tight- 
binding limit by a Bose-Hubbard model. The Hamiltonian of our system is given by 

H = J(r,r>t(r)a(r') + ^ ^ a^(r)2a(r)2, (1) 


where r and r' run over all the site positions of the lattice, a(r) is the bosonic 
annihilation operator at the site r, and U describes the on-site Hubbard interaction. 
The diagonal element J(r, r) gives a potential energy at the site r, and the off-diagonal 
element J{r,r') with r ^ r' describes the (generally complex) hopping amplitude 
between the two sites and satishes J{r',r) = J*{r,r') in order for H to be hermitian. 
We set the total number of particles to N: 

^a^(r)a(r) = W (2) 


We focus on the case in which the hopping terms in Eq. (|^), which will be denoted 
by hfkin, are given by the Haldane model in a honeycomb lattice |^; see Fig. |^. A 
honeycomb lattice consists of A and B sublattices, and is spanned by the primitive 
vectors ai ,2 («-3 is also introduced for convenience). We also introduce the vectors ^ 1 , 2 , 3 , 
which are directed from a B site to the three neighboring A sites. These vectors are 
given by 


0.1 = 62 — S3, 02 = S3 — ( 5 i, 03 = 5 i — S2, 


( 3 ) 
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Figure 1. (a) Haldane model on a honeycomb lattice. The kinetic part of the model, 
^^kin, consists of the nearest-neighbor hopping Ji, the next-nearest-neighbor hopping 
J 2 , and the potential difference 2A between the two sublattices A and B. Furthermore, 
atoms acquire an Aharanov-Bohm phase $ when hopping along every dashed line in 
the arrowed direction [see Eq. (||)]. (b) The first Brillouin zone of the honeycomb 
lattice. Here a\ and represent the reciprocal lattice vectors [see Eq. i)]- 


where d is the length between the neighboring sites. We also introduce the reciprocal 
lattice vectors 


= 


2n 

~d 


1 1 
3’71 


* 

an = 


271 

~d 


- -) 
3’73;’ 


(4) 


which satisfy ai ■ a* = 27i6ij {i,j = 1, 2). 

The Haldane model consists of the real nearest-neighbor hopping Ji, the complex 
next-nearest-neighbor hopping and the potential difference 2A between the two 

sublattices. Nonzero values of J{r,r') are thus given by 


— J(r -b 6j, r) = —J{r, r -|- dj) = —Ji for r E B; 

— J{r + aj, r) = —J*{r, r + aj) = — for r G A = A, B] (5) 

— J(r, r) = exA + V{r) for r G A = A, B, 

where j = 1, 2, 3, and eA,B = ±1. Here, r G A indicates that the site r belongs to the 
A sublattice, and V (r) describes an external potential which depends on the setting of 
our system. We assume Ji, J2 > 0 in the following. 


2.2. Band structure in the noninteracting case 


Here we set U = 0, and review the band structure of the Haldane model in 
noninteracting case 0 . Assuming the periodic boundary conditions in the 


directions of the honeycomb lattice, we perform the Fourier expansion 


the 

two 


a(r) 





(r G A 


A,B), 


( 6 ) 


where the sum is taken over the discrete momenta k in the hrst Brillouin zone, and Au^ 
is the total number of unit cells in the system (i.e., half of the total number of sites). 
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Figure 2. (color online) (a) An example of a band structure in the noninteracting case, 
calculated along the path Mi —>■ K- —>■ Kj^ —^ Mi in the first Brillouin zone shown in 
Fig. ^b). The parameters are chosen to be in the trivial phase with C± = 0. (b, c) 
Examples of Bogoliubov excitation bands in the presence of interaction. A transition 
in the band topology occurs at Un/ Ji = {Un/ Ji)c — 0.321 [a solution to Eq. (|^)], 
at which the band gap closes at the K- point [see (b)]. For Un/Ji > (t7n/Ji)c, the 
higher band acquires a nontrivial Chern number (7+ = +1. 


The kinetic part of the Hamiltonian (|l]) is then rewritten as 

= 'E (“h*:). A(fe)) nk) . (7) 

Here the 2x2 hermitian matrix 'H(fc) can be written in the form 

TLik) = ho{k)I + h{k) ■ cr, (8) 

where I is the identity matrix, cr = (ai, <72, cxs) are the Pauli matrices. The coefficients 
ho^k) and h{k) = {hi{k), h 2 {k), h^^k)) are calculated as 

ho{k) = — 2 J 2 cos(<h) cos(fc ■ ttj), hi{k) = —Ji cos(fc ■ 6j), 

' ' (9) 

h 2 {k) = —Ji sin(A; ■ <5^), h^i^k) = A + 2J 2 sin(<h) sin(A; ■ aj). 

3 3 

The two energy bands are obtained through the diagonalization of Eq. as 

e±(fe) = ho(fc) ± h(A;), h{k) := \h{k)\. (10) 

An example of energy bands is presented in Fig. @(a). For noninteracting fermions, 
complete filling of the lower band leads to a band insulator. 

For noninteracting bosons, Bose-Einstein condensation into the lowest-energy 
single-particle state occurs at zero temperature. || When J 2 = A = 0, the bottom of 


I In the thermodynamic limit, Bose-Einstein condensation does not occur at finite temperatures in 
two dimensions. In finite-size systems, however, a large condensate fraction can still be achieved if 
coherence is formed over the system at sufficiently low temperatures. 
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the lower band e_(fc) is located at fc = 0. To find whether the position of the bottom 
can change owing to finite J 2 or A, we expand Eq. (p^ aronnd fc = 0 as 


9 

e±(fc) = ^o(O) + cos <h ± 



9 \ 

4 ^9 + A2 J 


+ 0{k^). 


( 11 ) 


Therefore, e-{k) is minimized at fc = 0 if the following condition is met: 

J2 

— . ^ = + 2 J 2 cos d> > 0 . ( 12 ) 

This condition is satisfied in the regime J 2 , |A| <C Ji, which is relevant to the Haldane 
model realized in the scheme of Ref. [pO|] .P 

To determine the single-particle gronnd state, we parametrize h{0) nsing the polar 
coordinates as 


h{0) = (—3Ji,0, A) = h(0)(sin 6^0 cos TT, sin 6^0 sin TT, cos 6 * 0 ). (13) 

The k = 0 part of hfkin is then diagonalized by means of the transformation 

where we have defined the 2 x 2 unitary matrix 

JJ(Q -^^iI+a^)|2 -^9a,|2 ^ COs(0/2) sin(0/2)\ 

^ sin( 0 / 2 ) cos( 0 / 2 ) )' ^ ^ 

For noninteracting bosons, Bose-Einstein condensation occurs in the mode created by 
al(0). For interacting bosons, the condensate wave function is gradually modified with 
increasing the interaction, as discussed in the next section. 


3. Bogoliubov theory and excitation band topology for homogeneous 
condensates 


In this section, we present the Bogoliubov theory [^, for homogeneous condensates 
with weak repulsive interactions U > 0, and determine the band structure of Bogoliubov 
excitations. Here by “homogeneous”, we refer to the situation in which the system has 
the periodicity of the honeycomb lattice (we do not require the equivalence of the two 
sublattices). We then analyze the topology of the Bogoliubov excitation bands, and 
determine the parameter ranges that exhibit nontrivial topology. 


§ When A = 0 and $ = tt, the condition in Eq. ( 0 ) is written simply as Ji > 6 J 2 . By relating phases 
of bosons to angles of classical XY spins, we find that this condition is equivalent to the known stability 
condition of the ferromagnetic order in a honeycomb lattice magnet with competing ferromagnetic 2Ji 
and antiferromagnetic — 2 J 2 couplings [^, Q. 
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3.1. Bogoliubov theory 

To formulate the Bogoliubov theory for the present system, we first need to determine a 
condensate wave function by using the Gross-Pitaevskii (GP) theory. In the GP theory, 
we introduce the GP energy functional E by replacing {a{r),a^{r)) by (■^(r),'^*(r)) in 
the Hamiltonian (|lD, and minimize it with respect to {gp{r),'ip*{r)) under the constraint 
Ylr IV’(’’)P = Since the single-particle ground state is formed at fc = 0 (as discussed 
in Sec. we introduce the following homogeneous ansatz for the interacting case: 


^p{r) = 


VK. 


{TEX). 


(16) 


We also introduce the chemical potential /i as a Lagrangian multiplier to satisfy the 
particle-number constraint. The functional to be minimized is then given by 


E - ^,N ^ - ^,I] 



( 17 ) 


Minimizing this with respect to 'ipx -^) gives a homogeneous version of the GP 

equations: 


|W(0) 



= 0 . 


(18) 


Since the single-particle ground state is created by aL(0) in Eq. (|^, it is convenient to 
parametrize ('0 a,'0b)^ as 




Vn 


sin(0/2)\ 
cos{6/2)J ’ 


(19) 


where 0 = do when U = 0. Multiplying Eq. (|T^ by (/a, Jb) or (—/s, /a) from the left, 
we obtain 


ho(0) - h(0) cos{9 - 9o) + 2Un{f\ -F f^) = /r, (20a) 

h{0) sm{9 - 9o) + 2UnfAfB{fA - /I) = 0, (20b) 


where n := iV/(2iVuc) is the average number of particles per site. The parameter 9 is 
determined by solving Eq. (|20b|) ; the chemical potential fi is determined by substituting 
the obtained 9 into the LHS of (|20a|) . For the obtained 9, we consider the unitary 
transformation 


0 ^( 0 ) 

0 ^( 0 ) 


U{9,n) 



-Ib fa+\ 

fA Jb) \a-) 


( 21 ) 


Bose-Einstein condensation occurs in the lower-band mode created by aL. When A = 0, 
the two sublattices are equivalent, and thus Eq. ( pOb|) gives 9 = n/2 For |A| <C Ji, 
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we can expand Eq. (pOD in terms oi 9 — tt/2, obtaining 

A 


0 = ^ 
2 


+ o (aVj') , 


3 Ji + Un 

p = —3 Ji + Un — 6 J 2 cos ^ + O (A^/J) . 


(22a) 

(22b) 


As seen in Eq. (|22a|) , the potential difference 2A induces a density imbalance between 
the two sublattices; this imbalance is reduced by a repulsive interaction U > 0. 

We now discuss excitations from the condensate ground state by using the 
Bogoliubov theory. To this end, using Eqs. (|) and (pT]), we decompose a(r) into the 
condensate and noncondensate parts as 


air) = 


fx 

VK 


-.a- + a(r) (r E X). 


(23) 


Here, the noncondensate part is given by 

1 " 


a[r = 




^xfx(^+ + 'y ^ o-x{k)U'^ ^ 


{TEX) 


(24) 


with A = B and B = A. Following the Bogoliubov approximation, we replace both a_ 
and al by y/N, substitute Eq. (p^) into H — fiN, and expand H — up to quadratic 
order in a(r). The terms linear in a(r) or af(r) disappear because of the stability 
condition of the condensate in Eq. (|20b|) , and we arrive at the Bogoliubov Hamiltonian 


H — jjiN = - cA{k)M.{k)c)L{k) H— (a\_, a+'j ) “*" const. 


with 


ot\k) := (^a\{k),a^^{k),aA{-k),aB{-k)j . 

Here, we have introduced the 4x4 matrix Ai{k) and the 2x2 matrix A4+ as 


nik) - pJ + AUnF"^ 2UnF‘^ 

2UnF^ n^{-k) - ial + AUnF‘^ j ’ 


M{k) = 

M+ = [ho(0) + h(0) cos( 6 ' - 6 ^ 0 ) - p + Wnflfl] ^ + ^Unflf^ai 


(25) 

(26) 

(27) 

(28) 


with F := diag(/A,/B). 

To diagonalize the Bogoliubov Hamiltonian (| 
Bogoliubov transformations 


we perform generalized 


(x{k) = W{k)l3{k), 


a+ 


= 


K( 0 ) 

,h( 0 ) 


with 


/3t(fe) := (bl(k),bl(k),b4-k),b-(-k) 


(29) 


(30) 














CONTENTS 


10 



Figure 3. The Chern number C+ of the higher excitation band in the Haldane- 
Bose-Hubbard model with J 2 , |A| <C Ji- The boundaries between regions of different 
C+ values are given by Eq. (^, and indicated for Un/Ji = 0,0.5,1. The solid 
curves correspond to the noninteracting case . The regions with nontrivial topology 
(C+ ^ 0) enlarge as Un/ Ji increases. 


Here, W{k) and W+ are paraunitary matrices which satisfy 

W\k)T3W{k) = W{k)T3W\k) = rg, Wla3W+ = W+asWl = ag (31) 

with Ta := diag(l, 1, —1, —1). These equations ensure the invariance of the bosonic 
commutation relations. If the matrices W{k) and IV+ are chosen to satisfy 

WHk)M(k)W(k) ^ dmg(E+{k),E.{k),E4-k),E_{-k)), 

IVlM+W+ = E+(0)I, 

the Bogoliubov Hamiltonian (pSD is diagonalized as 

H-fiN = {k)bl(k)b^.{k) + (fc)6L(fc)6_(fc) + const.. (33) 

k k^O 

The paraunitary matrix W{k) satisfying Eq. ( [5^ ) can be constructed numerically by 
using the method described in Refs. Examples of the calculated Bogoliubov 

excitation bands E±{k) are presented in Fig. ||(b,c). 


3.2. Band gap 

Before discussing the topology of the Bogoliubov excitation bands, let us analyze the 
gap between the two bands E±{k). The band topology cannot change as far as the band 
gap remains open. Thus the closing of the gap can signal a change in topology. 

As seen in Fig. § and known in the noninteracting case the smallest gap is 
found at one of the two points k = EK := ± (a* + a^) /3 in the Brillouin zone; see the 
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K± points in Fig. ^ At these points, because of hi{EK) = h 2 {EK) = 0, the matrix 
Ai{-EK) in Eq. (p7|) is decoupled into A and B sublattice blocks, each of which can be 
diagonalized by a standard 2x2 Bogoliubov transformation. The excitation energies at 
k = EK are then calculated to be 


Ex(±lT) = Ax ± 3\/3ex J 2 sin $ (X = A,i?), (34) 

Ax := [(3 J 2 cos <!> + ex A - + AUnflf - {2Unflf] . (35) 

The higher (lower) of these energies gives E+(±1T) (E_(±1T)) in Eq. (^). The gap¬ 
closing conditions at fc = EK are thus obtained as 

0 = Ea{±K) — Eb{±K) = Xa — Xb A: 6\/3J2 sin 4). (36) 


For J 2 , |A| <C Ji, by using (^2]) , we can expand Xa — Xb as 


Xa — Xb — 


2A 


G{Un/Ji 


+ 0 


Jj ^ 
Ji Ji 


, G{s) ( 1 -|- 


2s 


1 + 


Equation 


is then rewritten into a simple form 

A ± 3 V 3 J 2 G{Un/Ji) sin $ + C> AV-^i) = 0. 


(37) 


(38) 


For Un = 0, this reduces to the exact phase boundaries A = =f3 \/3J2 sin <F in the 

indicates that with increasing the strength of 


noninteracting case Equation 

interaction Un/ Ji, these boundaries are shifted to larger |A| by a factor of G{Un/Ji)] 
see Fig. ^ Here, G{s) is an increasing function of s, and expanded for |s| 1 as 


G(s) = 1 + 2s/3 + sVl8 + 0{s‘^). 


(39) 


Across the gap-closing lines, the topology of the Bogoliubov excitation bands changes, 
as we discuss next. 


3.3. Chern number 

We now analyze the topology of the Bogoliubov excitation bands. We hrst note that 
technically, M(k) in Eq. (^Tl) (more specihcally, 'H(fc) in it) needs a modihcation 
for such a purpose because it is not periodic in the hrst Brillouin zone in the 
present representation. This is because the Fourier expansion (P) was based on the 
real-space positions r — while this treatment was useful in making ?{(k) possess 
the G 3 symmetry of the original lattice, the spacing between the two sublattices 
introduced an additional phase factor, which broke the periodicity in the hrst Brillouin 
zone. To recover the periodicity, we dehne A4{k) by replacing 'H(fc) by TL^k) = 
g-*fc-53(/-(T3)/2-^^^^gifc-(53(7-o-3)/2 (pTI). We then introduce the paraunitary matrix 

W{k) as the matrix which “diagonalizes” Ai{k) in the sense of Eq. (^2]). 
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The 4x4 paraunitary matrix W (fc) can be parametrized as 

^U{k) V*{-k) 


W{k) = 


V{k) W{-k) 


with 


lJ(h\ - (^A-{k)\ _ fVA+{k) VA-{k)\ 

^ ^ [uB+ik) UB-{k))^ ^ \VB+{k) VB-{k))- 

To discuss the topology of each excitation band, we introduce the vectors 

\w-,ik)) = {uA-,{k),UB-,(k),v^-,{k),VB;(k))'^ (7 = ±), 


(40) 


(41) 


(42) 


where 7 = + and — correspond respectively to the hrst and second columns of W{k). 
It follows from Eqs. ([TID and (^) that these vectors satisfy the eigen equation 


M[k)\w^{k)) = E^{k)T 3 \w-y{k)) 
and the orthonormality condition 


{w^{k)\T3\wY{k)) = 6^y. 


For each band, we introduce the Berry curvature [p7|, ^ 


= ieij{diW^{k)\T3\djW^{k)) 


(43) 


(44) 


(45) 


with di := The Chern number can then be introduced as [g, ^ 


= 


'BZ 


(Pk 


B^{k). 


(46) 


In the noninteracting case Un = 0, both C± are quantized to integers, and satisfy the 
zero sum rule — 0- interacting case Un > 0, however, \w_{k)) is not 

dehned at fc = 0,[[]] and thus there is an ambiguity in the dehnition of C_. Nevertheless, 
C+ is still well-dehned, and quantized to an integer.p| We numerically calculate (7+ 
using the manifestly gauge-invariant method proposed in Ref. [Q. 

The “phase diagram” based on the Chern number (7+ is presented in Fig. ^ (we note 
that this diagram is not based on ground-state transitions). We numerically conhrmed 
that the boundaries between regions of different C+ values are given precisely by the gap¬ 
closing condition (|^ (or Fq. (^^ for J 2 , |A| <C Ji) obtained in Sec. The obtained 


II This is because the Bogoliubov excitations consist only of modes orthogonal to the GP ground state 
|||§;see Eq. 

^ This does not contradict the zero sum rule. In the interacting case Un > 0, the rule applies to the 
sum over all the particle and hole bands. Namely, C+ + C- + -I- C'_ = 0, where C' 3 . are the Chern 

numbers associated with the hole bands . One can easily show C+ + C'j^ = 0 = C- + C '_, and thus 
the sum rule is trivially satisfied. Therefore, one cannot use the ill-defined nature of C- to change C+ 
to an arbitrary value. 
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results indicate that the topology of the Bloch bands known for the noninteracting 
case f/n = 0 are smoothly carried over to that of the Bogoliubov excitation bands 
for the interacting case Un > 0, and that the regions displaying nontrivial topology 
7 ^ 0 enlarge with increasing Un/Ji. When (7+ 7 ^ 0, the bulk-edge correspondence 
0 dictates that in-gap chiral edge modes intervening between the upper and lower 
bulk bands appear when the system has a boundary. We numerically demonstrate the 
emergence of such modes in Sec. |^. 

In closing this section, two remarks are in order. 

The hrst remark is on the reason why the ranges displaying topological bands 
expand with increasing the interaction U. In the noninteracting case U = 0, the 
potential difference 2 A between the two sublattices drives a transition from topological 
to trivial bands as it favors sublattice-separated Bloch wave functions, which have trivial 
topology. Algebraically, this potential difference induces a hnite difference between 
and Xb dehned in Eq. (^), and the transition occurs when |Aa — As| = 6 \/ 3 J 2 | sin<I>|. 
The interaction U > 0 has the effect of obscuring this difference [through Unfx in 
Eq. (^)], and thus a larger |A| is required to drive the topological-to-trivial transition. 
It will be interesting to investigate whether a similar stabilization of topological bands 
occurs in a wider variety of interacting systems. 

The second remark is on the case of attractive interactions U < 0. While an 
attractive Bose gas is unstable against collapse in the thermodynamic limit, it can form 
a metastable condensate in a hnite system if N is below a certain critical value . As 
far as a quasi-homogeneous condensate is realized, we can perform the same Bogoliubov 
analysis as in this section, and obtain the phase boundaries of topologically nontrivial 
regions as in Eq. (|38D; these regions gradually shrink with increasing \U\n for U < 0. 
In these regions, in-gap chiral edge modes discussed in Sec. ^ are also expected to be 
formed. 


4. Ground state and excitations in trapped condensates: formalism 


In this section, we present the formalism for calculating the ground state and excitations 
of trapped condensates. We hrst describe the Bogoliubov-de Gennes (BdG) theory 
1^ , ^ for inhomogneous condensates on lattices. We then apply this theory to a strip 
geometry, which is convenient for discussing edge modes. We also describe an extended 
Thomas-Fermi approximation which can give a simple analytic expression for the density 
prohle in a given trap potential. Numerical results obtained using the formalism are 
presented in Sec. ||. 


4 . 1 . Bogoliubov-de Gennes theory for inhomogeneous condensates 

The BdG theory for inhomogeneous condensates can be derived by linearizing a time- 
dependent GP equation. We start from the Heisenberg equation of motion for the 
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time-dependent operator a{r,t)- 

ihdtair, t) = [a(r, t),H] = — ^ J(r, r')a{r', t) + Ua^{r, t)a{r, t)‘^. (47) 

r' 

Replacing by classical fields ('^(r,t),t)), we obtain the GP 

equation 

ihdtijir, t) = — ^2 t) + U'ilj*{r, t)ip{r, t)^. (48) 

r' 

Equation (|^) imposes the normalization condition X]rlV’(’"G)P = Inserting a 

stationary state ansatz ipirO) = into Eq. (^Sj), we obtain the time- 

independent GP equation 

Ti^{r) = - E J{r,r')'ifj{r')+ Ui/j*{r)i/j{r)‘^. (49) 

r' 

The solution 'ipis') with the lowest frequency n/h gives the GP ground state. 

We now discuss small fluctuations around the GP ground state: 'ipirO) = 
'0(r)e“*^*/^ -f (f){rO)- Expanding the GP equation (^ 8 ]) to first order in 0(r), we obtain 
a linear differential equation 

ihdt(j){rO) = - ^ J(r,r')0(r',t) + 2f/|^/>(r)p0(r, t) + V(’^G)- ( 50 ) 


„iEt 


Assuming a solution of the form 

4>{r,t) = [u{r)e~^^^^^ + v*[r)e'^ 

we obtain Bogoliubov-de Gennes (BdG) equations 

+Eu{r) = — ^ J(r, r')u{r') + {2U\i/j{r)\‘^ — ii)u{r) + Ui/j{ryv{r), 

r' 

—Ev{r) = — ^ J*{r, r')v{r') + ( 2 f/|' 0 (r)p — iT)v{r) + U%1)*{rYu^r). 


(51) 


(52a) 

(52b) 


If the condensate is stable, these equations admit Ng — 1 (linearly independent) sets of 
solutions {uj{r),Vj{r)) with positive frequencies Ej/h, where Ng is the total number of 
lattice sites. Such solutions can be chosen to satisfy the orthonormality condition 

J2[u*jir)ui{r) - v*{r)vi{r)] = 5ji. (53) 

r 

The present formulation of the BdG theory is based on the linearization of the GP 
equation (48), which is an equation of motion for the classical field. However, we can 
check the consistency of this classical-field formulation with the operator formulation 
for the homogeneous case in Sec. |[ Indeed, substituting the ansatz (pTl) into the GP 


(r e X) 


(54) 


equation (^9|) reproduces Eq. (|T8|) . Furthermore, substituting 

•MZ 

into the BdG equation (^2D reproduces Eq. (^3]). It is known that the operator 
formulation for the inhomogeneous case also leads to the same set of BdG equations 
as in Eq. 


55 
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Figure 4. Strip (or “nanoribon”) geometry with zigzag edges. A periodic boundary 
condition is imposed only along the y direction. The red square indicates a unit cell 
which contains = 16 lattice sites. We denote the number of unit cells by Ny. The 
origin of the x coordinate is placed at the center of the system. 


4-2. Strip geometry 


We apply the BdG theory to analyze excitations in trapped condensates. We describe 
site positions by two-dimensional coordinates r = {x,y). For simplicity, we consider a 
strip geometry, which is periodic only along the y direction as in Fig. ^ and analogous 
to a graphene nanoribbon. This geometry can be used to describe the central part 
of an elongated condensate around which the system is approximately uniform in the 
elongated direction. Along the x direction, we introduce a box trap with sharp zigzag 
edges (Fig. |) or a harmonic trap in Sec. however, we do not assume a specific 
trap potential in this subsection. 

Exploiting the translation invariance along the y direction, we make the following 
ansatz for the GP ground state: 


V’(x,?/) = 






y/'^yfxi 


:= 


N 


(55) 


where Ny is the number of unit cells, and i/j^ = \/N satishes the normalization 
condition Xlx iV'xP = ~ Substituting this into the GP equation (^91) , 

we obtain 

^-'(0)/-' + Uriyrjl (56) 


Here we have introduced 


y-xx' (ky) = a;', y') 


,-iky{y-y') 


(67) 


where {x', y') is a particular site position, and the sum over y is restricted to the site 
positions for hxed x; the translation invariance along the y direction ensures that the 
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RHS of Eq. does not depend on y'. An accurate solution to Eq. ( ^61) with the lowest 

frequency y/h can be obtained by numerically performing the imaginary time evolution 
with the replacement y —dr- 

After obtaining the GP ground state (|^, we solve the BdG equations to 
calculate excitations. We introduce the following ansatz with momentum ky in the y 
direction: 


u{x,y) = 


He 


v{x,y) = 


Jkyy 


Substituting these into the BdG equation (|5^, we obtain an eigen equation 


(58) 


M.[ky)w = Et^w. 


(59) 


Here we have introduced a 2A3;-component vector w = {{ux}, and a 2Nx x 2Nx 

matrix 


Miky) = 


Tiiky) — yl + 2UnyF*F 


UUyF^ 


UuyF*^ 


-tI + 2UnyF*F / ’ 


(60) 


with F := diagd/a;}). The eigen equation (|59|) can be solved numerically by 
using the method of Ref. |^. We denote positive-frequency solutions by Wj{ky) = 
{{uxjiky)}, {vxjiky)})^ and Ej{ky) with j = 1, 2 ,..., A^. in descending order in energy. 
In the operator formulation of the BdG theory, this leads to the fact that the 
Hamiltonian is diagonalized as 

H — fiN = Ej{ky)b^j{ky)bj{ky) + const., (61) 

(ky,j)rtiO,Nr) 

where the new bosonic operators {bj{ky)} are related to the original ones as 

1 


a{x,y) = i>{x,y) + 




{kyJ)^{0,Nr) 


ikyV 


Uxjiky)bj{ky) + ky)bj( ky) . (62) 


To discuss excitations in a strip geometry, it is useful to introduce the spectral 
weight at zero temperature 


p{x,x',ky,u) = / dt'^e ^’^y^y y'^"^'''^Wa{x,y,t),a\x',y',C)]), 


(63) 


where {x,y) and {x',y') are taken similarly as in Eq. (^71), and the average (•) is taken 
over the vacuum of the Bogoliubov excitations in Eq. (|6TD. Using Eqs. (^ and 
Eq. (|63D can be rewritten as 


p{x,x',ky,u) = 27i'^[uxj{ky)ul,j{ky)6{u - Ej{ky)/h) 
j 

+ (-y)<Vj(-y)5(w + Ej(-ky)/K)]. 


(64) 


In Sec. we calculate p{x,x,ky,oj) to discuss excitations that can be probed at each 
position X. We note that similar calculations are performed for a fermionic Hofstadter 


model in various traps by Buchhould et al. 










CONTENTS 


17 


Ji^.3. Extended Thomas-Eermi approximation for density profiles 


Here we provide an approximate solution to the scaled GP equation (p6[). This helps us 
tune the potential and the interaction to obtain the desired density prohles later. We 
extend the Thomas-Fermi approximation ^ , where the condensate wave function is 
assumed to vary slowly in space, so that the oscillations of fx between the two sublattices 
due to A 7 ^ 0 are also taken into account. We assume 


fx = f{x) + ex5f{x), (65) 

where f{x) and 6 f{x) are slowly varying real functions, and ex = +1 (—1) if x belongs 
to the A (B) sublattice. As can be seen in Fig. all the sites having the same x belong 
to the same sublattice. Assuming |A| A, we can expect \6f{x)\ <C f{x). Then the 
RHS of Eq. (|56D can be approximated as 

a:+3/2 
x'=x—3l2 

^ [V(x) — 3 Ji — 6 J 2 cos $ + Unyf{x)'^] f (x) 

+ ex [A/(a;) + (V{x) + 3A - 6 J 2 Cos<l> + 3Unyf{x)‘^)6f{x)] , 

where we have ignored higher-order terms in A and df{x). Here, the hrst and second 
lines of the last expression can be viewed as uniform and staggered components because 
ex oscillates rapidly and other functions of x vary slowly. Requiring that Eq. (|56D holds 
separately for different components (because they cannot cancel each other), we obtain 


f{x) A 

a;+3/2 


x'=x—3l2 


6 f{x) + Uny[f{x) + 6 f{x)f 


r66'l 


Unyf{xY = max(t/n““ — V (x), 0), := /u -h 3 Ji + 6J2 cos $, 


df(x) = 


Af(x) 


V(x) — fi -h 3Ji — 6 J 2 cos <F -|- 3Unyf{xy 
Therefore, the density prohle scaled by the interaction U is obtained as 

Unylfxl"^ ~ max(7/nT^ — Vix)^ 0) ( 1- . - 

^ -V{x) + 3Ji 


(67a) 

(67b) 


( 68 ) 


This expression agrees well with the density prohles of the numerically calculated GP 
ground states presented in Sec. |^. As seen in this expression, can be interpreted 
as the maximal uniform-component density achieved at the potential minimum with 
V{x) = 0. 


5. Ground state and excitations in trapped condensates: numerical results 

In this section, we present numerical results on the ground state and excitations of 
trapped condensates based on the formalism described in Sec. We exploit the strip 
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Figure 5. Density profile (scaled by U/ Ji) of the ground state along the x direction 
in a box trap with = 80 (see Figjl). The data for the GP ground state are 
compared with the extended TF result (^. The model parameters are chosen such 
that the case of Fig. |2|(c) with interaction-induced nontrivial topology = -1-1 is 
realized in a quasi-homogeneous region in the bulk. 


geometry described in Sec. [4.2|, and introduce a box trap or a harmonic trap along the 
X direction. While harmonic traps are used more commonly in experiments of ultracold 
atoms, a box trap with sharp boundaries has also been realized recently |^. Sharp 


boundaries can also be realized in synthetic dimensions [26 


], We demonstrate 
that chiral edge states reflecting the nontrivial bulk band topology do appear in a box 
trap. For a harmonic trap, by contrast, our results show that such edge states are 
substantially obscured and difficult to observe, as opposed to an expectation from a 
semiclassical picture. 


5.1. Box trap 


We first consider the case of a box trap 


V{x) 


0 (|a;| < \N^d) ; 

oo (otherwise). 


(69) 


which has sharp boundaries of zigzag type as in Fig. ^ The extended TF result (^) 
leads to a uniform density in each sublattice inside the box. Summing Eq. (|68D over the 
lattice sites |x| < ^Nxd, we find lluy ~ Urd^^Nx. We note that in the BdG calculation 
for a strip geometry, the product Uriy is an input parameter, i.e., if llUy is fixed, the 
result does not depend on individual values of U and Uy. We can thus tune Uriy to 
obtain a desired scaled average density Urd^^/Ji ~ llUy/(N^Ji). 

We set {J 2 e^^/ Ji, A/ J 2 ,Uny/(N^Ji)) = (O.li, 1.2 x 3^/^, 1) so that the case of 
Fig. |^(c) is realized in a quasi-homogeneous region in the bulk. The scaled density 
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(a) X interval: (-30d,-15d) (b) x interval: (+15d,+30d) 



Figure 6. Integrated spectral weight (^) xini = (—30d, — ISd) 

and (b) a; in / = (+15(i,+30d), which contain the left and right edges, respectively. 
The BdG calculation is performed using the GP ground state shown in Fig. ||, and the 
spectral weight p{x,x,ky,uj) is calculated with Eq. (^. Positive (negative) energies 
correspond to particle (hole) excitations. 


profile Uriylfxl'^/Ji of the GP ground state in Fig. is indeed almost uniform in each 
sublattice, and agrees well with the extended TF result (|68D , except over a few sites 
near each boundary. 

We now discuss excitations calculated by the BdG theory described in Sec |^ . 
For an interval I of the x coordinate, we consider the integrated spectral weight 
'n,x£i 1^)1 where p{x,x,ky,u) is dehned in Eq. (|63|). The results for (a) 

/ = (—30d, —ISd) and (b) I = (+15d,+30(i) are presented in Fig. In both of 
these results, the continuum of excitations corresponding to the two bulk particle (hole) 
excitation bands is found with high (low) spectral density. Furthermore, inside the band 
gaps, chiral modes with negative and positive velocities are clearly formed in (a) and 
(b), respectively.|3 These results are consistent with the formation of chiral edge modes 
propagating in the —y {+y) direction at the left (right) boundary as expected from the 
bulk topological number (7+ = +1. 

To detect the chiral edge modes between the two excitation bands, a high-frequency 
probe is required. This contrasts with the case of fermionic topological insulators, 
where edge modes cross the Fermi level and can be excited with inhnitesimal energies. 
In ultracold-atom experiments, stimulated Raman transitions can be used to create 
excitations with desired momentum and frequency. Since edge modes are isolated from 
bulk modes in momentum and frequency in Fig. p, Raman transitions can transfer a 


+ In Fig. ^ [and Fig. ||(c)], the bulk band gap is relatively small because we examine the case in which 
nontrivial topology G+ = -1-1 is induced by the interaction U. That is, the system is located in a narrow 
region between the lines for Unj J\ — 0 and 1 in Fig. where a small band gap closes and opens again 
as we change Unj J\ between these values. If the system is located deep inside a region with nontrivial 
topology, a much larger band gap can be created, and edge modes can then be distinguished more 
clearly from bulk modes. 
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Figure 7. Interference patterns of an edge matter wave with the background 
condensate at different times Jit/h = 0,0.5,1,1.5, 2. The color shows the scaled 
density relative to that of the ground state, [|V'(x, y, |/a,p], calculated 

from Eq. dnl)- The center of each triangular pixel corresponds to a site of the 
honeycomb lattice. The edge matter wave is created by transferring a portion of 
the condensate into the edge mode with the momentum ky{'/5d)/{2^) = —0.35 
(equivalent to 0.65 due to the periodicity of the Brillouin zone) and the energy 
E 4 o{ky)/Ji = 4.03 as can be seen from Fig. ||(b). During the time evolution (a)- 
(e), the pattern propagates in the (negative) y direction along the right edge with the 
phase velocity E 4 o{ky)/{kyJi) = —1.83 x v^d; see, e.g., the propagation of an antinode 
with a positive variation (red) as indicated by the arrows. 


portion of a condensate selectively into a particnlar edge mode with momentnm ky and 
freqnency u = Ej{ky)/h, realizing an “edge matter wave.”Q The resulting condensate 
wave function is a coherent superposition of the background condensate and the edge 
matter wave, which is given by 

y, t) ^ [f, + , (70) 

where a is the complex amplitude of the edge mode. We assume |ap <C 1 to ensure 
that the linearization done in Eq. ( 0 ) to derive the BdG theory works. We relate the 
amplitude a to the microscopic process later. Since the edge mode has its weight mainly 
around an edge, we can expect that a density wave is formed along the edge as a result 
of the interference with the background condensate. Using Eq. o. the scaled density 
prohle relative to the ground state is calculated as 

II Iln 

-r -^yl/xT] ~ [az^j{ky)E^^^^-‘^^^ + c.c.] (71) 

Ji Ji 

* While the edge matter wave has an infinite lifetime within the BdG theory, it can acquire a finite 
lifetime due to collisions between quasiparticles and condensed particles (known as Beliaev damping 
1^). The estimation of this lifetime however requires a detailed analysis of the collision processes, 
which is beyond the scope of the present paper. 
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with Zxj{ky) := f^Uxjiky) + fxVxj{ky). We shift the origin of t such that a becomes real. 
We plot Eq. for different times in Fig. |^. A density wave is indeed formed along 
the right edge, and it propagates in the negative y direction with the phase velocity 
EiQ{ky)/{kyJi) = —1.83 X \f^d of the edge mode. We expect that such a propagating 
density wave can be used as a macroscopically enhanced experimental signature of an 
edge mode. 

In the above argument, we have kept the amplitude a of the edge matter wave 
undetermined. Here we determine a based on the microscopic process. This would help 
design an experimental setup for creating and observing an edge matter wave. A pair of 
Raman lasers with wave vectors fci 2 and frequencies Ui ^2 ai'e prepared in such a manner 
that k = ki — k 2 points in the y direction [i.e., k = (0, ky, 0)], and that fuv = hiui — 002 ) 
is resonant with the excitation energy Ej{ky) of the edge mode. These lasers induce the 
following time-dependent perturbation to the Hamiltonian: 

H\t) = E hO,cos{kyy — ujt)a^ {x,y)a{x,y), (72) 

0 ,y) 

where the sum runs over positions {x, y) of lattice sites, and H describes the strength of 
the atom-light coupling. This perturbation adds the following term to the RHS of the 
linearized GP equation (|50[): 

Ml cos{kyy — ut)‘ip{x, y, t) ^ Ml cos{kyy — . (73) 

This describes a transfer of the condensate particles to the target edge mode; higher- 
order processes in which particles are further transferred to higher-energy states are 
neglected. In the presence of this term, we assume a solution of the form 

0(a;, y, t) = [a{t)ux,{ky)E^’^yy-^^^ + a*(f)nT(fc,)e-*('=-"--*)] . (74) 

The amplitude a{t) of the edge matter wave is found to obey 

[ihdtait)] T^Wj{ky) = (75) 

where Wj{ky) = {{uxj{ky)}, {vxj{ky)})'^ as dehned in Sec. and f := {{fx},{fx}V- 
If Raman lasers are illuminated in the time interval [0,hf], the integrated amplitude of 
the edge matter wave is obtained as 

a{6t) = -i^w]{ky)f, (76) 

where we have used wj{ky)T 3 Wj{ky) = 5ij in solving Eq. ([75|) . Since the GP ground 
state / is extended over the strip and the edge-mode wave function Wj{ky) is localized 
only over a few sites around the edge, their overlap scales as w^j{ky)f ~ l/\/Nx. Using 
this result, we can tune G or 6 t to achieve a desired amplitude a. 
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Figure 8 . (a) Density profile (scaled by 17/Ji) of the ground state in a harmonic 

trap. The data for the GP ground state are compared with the extended TF result 
(1^. The model parameters are chosen so that the case of Fig. ||(c) is realized around 
the center of the trap. Calculations are performed in a sufficiently wide strip (with 
Nj; = 160) so that the effect of the strip edges is negligible, (b) Local band edges 
Ea/b{—K) calculated from the density profile in (a). For each x G X, Ex{—K) is 
calculated by substituting the local density Uny\fx\'^ and the local potential V{x) — ^ 
into 2Unfx and —respectively, in Eq. (|4|); Ex{—K) is calculated with a similar 
procedure using the average density at the two neighboring sites. 


5.2. Harmonic trap 


We next consider the case of a harmonic trap V{x) = 

(|68|) suggests that the condensate extends over the range |a;| < i?TF 
Integrating Eq. (|68|) over this range, we hnd 


The extended TF result 


Uriy'GG 


dx{Rl^ - x^) = 


16 i?' 


■TF 


3d 


9d 


9d 




(77) 


—Rtf 

To achieve the desired values of and T/tf, we can thus set the input parameters 

as 

(78) 


2 


d2 ’ 

rurpp 


Urn, = 


9d 


'-TF 


We set j di, A/J2, / Ji) = (O.li, 1.2 x 3^/^, 1) so that the case of Fig. Kc) 

(with nontrivial topology C+ = +1) is realized around the center of the trap. The 
scaled density prohle Uny\fx\^jJi of the GP ground state in Fig. |(a) indeed shows 
the maximum of near unity at the center; the oscillating pattern of the prohle agrees 
well with the extended TF result (^). As we move away from the center, the density 
decreases towards zero, and the case of Fig. ^(a) (with trivial topology (7+ = 0) is 
expected to be realized for |x| > i?TF- This indicates that topological boundaries appear 
inside the condensate, assuming local homogeneity as in the semiclassical approach. 
To locate such boundaries, we plot in Fig. |^(b) the local band edges Ex{—K) with 
X = A, B, which are calculated by substituting the density prohle of Fig. ||(a) into 
Eq. (^). The two energies Ea/b{—K) indeed crosses at x/d ~ —26 and 24.5, where 
topological boundaries are expected to be formed. 
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Figure 9. (a) Local spectral weight p{x,x,ijj) = p{x,x,ky,uj). A BdG 

calculation is performed using the GP ground state displayed in Fig. ||, and the spectral 
weight p{x,x,ky,uj) is calculated with Eq. (^3|). (b) Integrated local spectral weight 
p(a:, a;, w) for small intervals I = (0,3d), (24d, 27d), and (30d, 33d). 


To see whether edge modes of the topological origin appear at such boundaries, 
we plot in Fig. |](a) the local spectral weight p{x,x,u) = Yky 
formation of an excitation gap with the vanishing spectral weight is seen around the 
center x = 0. As we move away from the center, the gap gradually closes as expected 
from the semiclassical approach. However, the reopening of the gap as in Fig. H(b) is 
not seen in this hgure. This is more clearly seen in the integrated spectral weight for 
small intervals shown in Fig. ^(b). For I = (0,3d), the formation of a gap with the 
vanishing spectral weight can be seen. For I = (24(i, 27d) and (dOd, 33d), by contrast, 
the spectral weight has a tiny but non-vanishing value at the valley, indicating a gapless 
nature. This indicates the breakdown of the semiclassical approach. A possible reason 
for it is as follows. The local band edges calculated with the semiclassical approach 
in Fig. P(b) indicate that not only the size of the gap but also its location (in energy) 
changes as a function of the position x. Therefore, beyond the semiclassical picture, 
the energy gap for a particular position x is easily penetrated by the states in the same 
energy range in the surrounding region. The observation of edge modes in a harmonic 
trap thus remains a challenging issue. 

6. Summary and outlook 

We have studied the topological properties of Bogoliubov excitation bands in BECs in 
optical lattices on the basis of a Bose-Hubbard extension of the Haldane model. We have 
shown that the topological properties of the Bloch bands in the noninteracting case are 
smoothly carried over to those of Bogoliubov excitation bands in the interacting case, 
and that the parameter ranges showing nontrivial topology enlarges with increasing the 
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Hubbard interaction or the particle density. In the presence of sharp boundaries, chiral 
edge modes appear in the gap between topologically nontrivial excitation bands. We 
propose that Raman transitions can be used to coherently transfer a portion of the 
condensate into an edge mode, and that a density wave is formed along the edge due to 
an interference with the background condensate. This can be used as a macroscopically 
enhanced experimental signature of the edge mode. By contrast, our results for a 
harmonic trap show that edge states are substantially obscured and difficult to observe, 
as opposed to what is expected from a semiclassical picture. 

We expect that BECs in optical lattices offer a unique playground in the studies of 
band topology. The macroscopic nature of BECs can enhance signatures of a topological 
edge mode. The high controllability of BECs offers various methods of exciting particles 
to such edge modes. While we have considered Raman transitions in this paper, a trap 
quench can provide another useful method. While both the bulk and edge modes are 
excited by such a quench, the edge excitations may exhibit a distinct time evolution 
because of their chiral nature (see Ref. for related numerical demonstrations for 
fermions). It will also be interesting to exploit the high controllability of optical lattices 
to design bosonic systems with different symmetries or dimensionality, where different 


topological classes can be explored as expected from the studies of fermions |^, |T8 
Note added. —Recently, we became aware of two independent works 0,0, where 
nontrivial topology of Bogoliubov excitation bands and associated edge states were 
discussed in different bosonic systems. Engelhardt and Brandes 0] have considered 
a one-dimensional system with inversion symmetry, while Bardyn et al. 0| have 
considered a kagome vortex lattice with potential realization in nonlinear optical systems 
or exciton-polariton condensates. 
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